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摘 Xm. 以 库 姆 塔 格 沙漠 南部 阿尔 金山 北 末 山 前 戈壁 区 为 研究 区 ,借助 无 人 机 影像 数据 准确 提取 植 
被 履 盖 区 ;采用 高 空间 分 辨 率 WorldView-3 数据 的 估算 结果 对 中 空间 分 辩 率 Landsat-OLI 数 据 的 全 
算 结 果 进 行 修正 ;选取 红 光 反射 率 波段 (RED)、 近 红外 反射 率 波段 (NIR)、 比 值 植被 指数 (RVI) VH 
一 化 植被 指数 (NDVI) 和 修正 土壤 调节 植被 指数 (MSAVI) 反 射 率 波段 作为 遥感 特征 参数 ,采用 逐步 
线性 回归 (SLR) 、 套 索 回 归 (LASSO) 和 岭 回 归 (RR) 模 型 分 别 对 稀疏 荒漠 植被 地 上 生物 量 进 行 估 
算 。 结 果 表 明 :(1) WorldView-3 数据 可 准确 刻画 极 干旱 区 稀疏 荒漠 植被 的 时 空 分 布 特征 ,其 估算 
精度 高 于 Landsat-OLI 数 据 。(2) 采用 非 线 性 的 多 元 回归 模型 SLR .LASSO 和 RR 构建 还 感 特征 参数 
与 地 面 样 方 数据 之 间 的 回归 模型 可 以 有 效 提高 模型 的 稳定 性 。RR 模型 的 回归 效果 优 于 SLR 模型 
和 LASSO 7 78 , (3) 提取 稀 政 荡 漠 植 被 信息 时 ,采用 WorldView-3 数 据 的 植被 生物 量 估 算 结 果 对 
Landsat-OLI 数 据 的 植被 地 上 生物 量 估算 结果 进行 修正 ,可 以 有 效 提高 Landsat-OLI 数 据 的 提取 精度 。 
X 键 词 : 遥感 特征 参数 ; 稀疏 荒漠 植被 ; 地 上 生物 量 ; 回归 模型 


文章 编号 : 


我 国 干旱 区 总 面积 约 为 2.15x10' km , 约 占 国土 
面积 的 22%。 包 括 新 疆 全 境 、 甘 肃 河西 走廊 地 区 和 
内 蒙古 阿拉 善 盟 及 宁夏 贺兰山 以 西 地 区 "。 匾 漠 植 
被 主要 分 布 于 干旱 区 ,以 旱 生 ` 超 旱 生 的 小 乔木 HE 
AR . 半 灌 木 和 小 半 灌 木 为 建 群 种 ,是 衡量 荒漠 生态 
系统 健康 状况 的 重要 指 征 一 。 

植被 地 上 生物 量 是 指 单位 面积 内 所 有 地 上 植 
物 的 总 重量 ,通常 以 鲜 重 或 干 重 表示 ,其 估算 方法 
一 直 是 生态 环境 监测 领域 的 研究 热点 ””。 传 统 的 
植被 地 上 生物 量 测算 主要 依靠 地 面 样 方 调查 ,虽然 
精度 高 ,但 费时 费力 .实施 困难 ,对 于 复杂 生态 系统 
存在 一 定局 限 ,只 适用 于 小 尺度 人 研究。 与 传统 方 
法 相 比 ,遥感 技术 在 植被 空间 异 质 性 识别 和 不 同时 
相 植 被 信息 获取 等 方面 的 优势 ,使 其 成 为 大 尺度 植 
被 地 上 生物 量 佑 算 的 主要 手段 。 


收 稿 日 期 : 2021-04-19; ”修订 日 期 : 2021-11-26 


近年 来 ,光学 遥感 植被 信息 提取 技术 被 广泛 
应 用 于 森林 和 草原 植被 ,但 是 在 荒漠 植被 中 应 用 
较 少 " 中 。 这 一 方面 是 由 于 蕊 漠 植 被 结构 简单 、 低 矮 
Fai .单位 面积 地 上 生物 量 非 常 低 ,传感器 对 稀 玻 
车 漠 植 被 的 敏感 性 弱 ,通用 遥感 模型 的 估算 精度 
低 ; 另 一 方面 在 提取 翘 漠 植被 信息 过 程 中 容易 出 现 
背景 (主要 是 土壤 ) 信 息 对 目标 信号 的 “污染 ”, 进 一 
步 降低 了 遥感 模型 的 估算 精度 。 为 了 解决 植被 信 
息 提取 精度 较 低 的 问题 ,Cunliffe 等 由 在 小 尺度 范围 
内 ,利用 无 人 机 遥感 影像 构建 改进 的 多 端 元 光谱 分 
解 模 型 来 提取 干旱 区 植被 覆盖 度 ,取得 了 较 好 的 效 
果 。Yan 等 "在 毛乌素 沙 地 采用 Landsat TM FIMO- 
DIS 遥感 数据 ,综合 分 析 了 归 一 化 植被 指数 (ND- 
VD 、 比 值 植被 指数 (RVI) 、 差 值 植被 指数 (DVI) 、 土 
二 调节 植被 指数 (SAVI) 修正 的 土壤 调节 植被 指数 
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(MSAVI) 5 种 植被 指数 与 植被 覆盖 度 的 相关 性 , 结 
果 显 示 MSAVI 与 稀 玻 的 沙 地 植被 覆盖 度 的 相关 性 
更 高 。 叶 静 共 等 ”采用 QuickBird 遥感 数据 ,基于 地 
面 调查 数据 ,对 乌 兰 布 和 沙漠 东北 缘 荒 漠 - 绿 洲 过 
渡 带 的 植被 地 上 生物 量 进行 了 估算 ,发 现 采 用 RVI 
线性 模型 估算 荒漠 植被 地 上 生物 量 效果 较 好 。 在 
植被 地 上 生物 量 模 型 方面 ,Zandler 等 1 利用 多 种 方 
法 对 塔吉克 斯 坦 帆 米尔 高 原 东部 稀 玻 荒漠 植被 的 
地 上 生物 量 估 算 结果 进行 比较 后 发 现 , 套 索 回 归 
(LASSO) 的 估算 效果 最 好 。 综 上 所 述 ,空间 分 辩 率 
更 高 的 遥感 数据 和 不 断 优 化 的 数据 分 析 方 法 可 以 
有 效 提高 稀疏 车 漠 植 被 遥感 信息 提取 的 精度 。 

为 了 提高 稀 玻 荒漠 植被 地 上 生物 量 的 估算 精 
BE ,本 研究 以 库 姆 塔 格 沙漠 南部 阿尔 金山 北 划 山 前 
戈壁 区 为 研究 区 ,从 遥感 数据 源 和 模型 方法 2 个 方面 
进行 探索 和 尝试 ,以 建立 更 加 准确 、 实 用 的 稀 玖 荒漠 
植被 地 上 生物 量 遥 感 佑 算 方 法 ,以 期 为 干旱 区 荒漠 
生态 系统 碳 储量 估算 和 荒漠 化 监测 提供 科学 依据 。 


1 研究 区 概况 


WEE PCT REA e LU JO EE REUS PE JS , 北 接 库 


姆 塔 格 沙漠 和 库 鲁 克 塔 格 山 , 总 面积 为 9397.25 km’, 
地 理 位 置 介 于 39°09'~39°50'N 与 90°55'~94°28 福 之 
间 ,海拔 1400~2900 m, 地 势 南 高 北 低 , 东 高 西 低 ( 图 
1)。 根 据 2017 一 2020 年 的 观测 ,研究 区 年 降水 量 大 
约 在 50~180 mm 之 间 , 在 空间 上 由 东 向 西 递 碱 ,由 
南 向 北 递减 。 降 水 多 集中 在 5 一 7 月 ,常常 形成 洪 
水 ,使 洪 积 房 上 沟 蜜 纵横 。 在 地 势 较 低 的 沟 床 内 主 
38 ^) fg 4 W E. (Zygophyllum xanthoxylon) WF XR 
(Gymnocarpos przewalskii ) MIF REE (Ephedra prze- 
walskii) 等 灌木 ,局 部 地 段 分 布 有 零星 的 盐 生 草 
(Halogeton glomeratus ) WIYE (Salsola ruthenica ) 和 
25 UK (Bassia dasyphylla) 等 一 年 生 草 本 植物 ,群落 
盖 度 大 多 在 10%~20% 之 间 ,灌木 植株 高 度 一 般 小 于 
1 m, 草 本 植物 高 度 一 般 不 超过 10 em ; 在 地 势 较 高 
的 沟 岸 地 段 主要 分 布 着 合 涉 获 (Sympegma regelii) 、 
红 砂 (Reaumuria soongarica) 等 灌木 ,群落 盖 度 在 
10% 以 下 ,许多 地 段 甚至 在 1% 以 下 ,植株 高 度 大 多 
为 10~50 cm; 在 洪 积 肩 北部 沙漠 边缘 地 带 分 布 有 梭 
T£ (Haloxylon ammodendron) , W $4 X (Calligonum 
mongolicum ) , PE Sl (Tamarix ramosissima ) 等 灌木 ， 
群落 盖 度 大 约 10%~20% ,植株 高 度 有 时 可 达 2.0~ 


图 例 

。 地 面 样 方 

O FAR 

CI WorldView-3 影 像 覆 盖 区 域 
O 无 人 机 影像 覆盖 区 域 

O WRK 


图 1 研究 区 位 置 以 及 样 方 分 布 
Fig. 1 Location ofthe study area and distribution of plots 
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2.5 mi"), HD ie PAE BAY FEL REY £8 JG VR 
多 样 性 和 生物 量 的 空间 异 质 性 均 处 于 较 高 水 平 ,对 
采用 遥感 方法 进行 植被 生物 量 的 准确 估算 提出 了 
挑战 。 


2 数据 与 方法 


2.1 地 面 样 方 调查 

2017 年 9 月 ,采用 样 带 与 样 方 相 结合 的 方法 进 
行 野外 调查 。 根 据 阿尔 金山 北 药 戈 壁 区 山 前 荒漠 
植被 分 布 的 垂直 地 带 性 特征 ,在 研究 区 的 东西 方向 
每 隔 20~50 km 设置 一 条 南北 方向 样 带 。 沿 着 样 带 
方向 每 隔 3 km 设置 一 个 30 mx30 m 的 样 方 。 调 查 
的 12 条 样 带 和 100 个 样 方 的 具体 分 布 (图 1)。 调 查 
内 容 包括 每 个 样 方 的 地 貌 特 征 .土壤 特征 ,海拔 和 
中 心 点 GPS 坐标 与 植被 生长 状况 ,测量 并 记录 样 方 
内 每 株 灌木 植物 的 基 径 (D,cm)、 冠 幅 长 轴 (C,cm) 
和 冠 幅 短 轴 (C;,cm) ARE (H, cm) ,以 及 每 个 植株 的 
相对 位 置 。 参 考 《 全 国生 态 状况 调查 评估 技术 规 
范 一 一 荒漠 生态 系统 野外 观测 》 中 的 藻 漠 植被 采 
样 方法 。 在 研究 区 内 分 不 同 植物 种 选取 若干 株 
不 同 大 小 的 健康 植株 ,分 别 测量 其 基 径 、 株 高 . 冠 幅 
长 轴 和 和 冠 幅 短 轴 ,然后 整 株 贴 地 剪 下 , 称 其 鲜 重 后 
取样 带 回 实验 室 ,在 80 C 通 风干 燥 箱 内 烘 干 48h 称 
重 , 经 计算 得 到 植株 地 上 生物 量 干 重 。 少 量 样 方 有 
草本 植物 分 布 ,在 有 草本 植物 分 布 的 样 方 内 ,在 样 
方 任意 1 条 对 角 线 的 两 端 和 中 间 各 取 1 个 1 mX1 m 
的 小 样 方 ,采用 收割 法 测定 小 样 方 内 所 有 草本 植物 


地 上 生物 量 鲜 重 , 并 采样 带 回 实验 室 , 在 80 通风 
干燥 箱 内 烘 干 48h 称 重 , 经 计算 得 到 每 个 小 样 方 草 
本 植物 地 上 生物 量 干 重 。 
2.2 样 方 地 上 生物 量 估 算 

异 速生 长 模型 法 是 估算 灌木 地 上 生物 量 的 有 
效 方法 ,可 以 减少 对 植被 的 破坏 ""。 基 于 地 面 样 方 
调查 数据 ,分 别 以 每 种 灌木 的 基 径 (D,cm) 、 冠 面积 
(A=1C,C,/4, cm’) AR (H, cm) FERRIERE (VAH, 
cm ) 为 预测 变量 ,构建 异 速生 长 方程 ,并 选取 决定 系 
数 ( 尼 ) 和 平均 绝对 误差 (MAE ) 作 为 比较 和 评判 异 速 
生长 方程 拟 合 效果 的 指标 。 计 算 公 式 如 下 : 


e=1-$o,-1 / 3-5 (1) 


MAE- 797-3: (2) 

式 中 :为 样本 个 数 ; y 为 单 株 灌木 的 生物 量 实测 值 
(g); 9, 为 单 株 灌木 的 生物 量 估 测 值 (g); 元 为 每 种 
灌木 的 生物 量 实测 平均 值 (g)。 

比较 估算 结果 ,每 种 灌木 的 异 速 生长 方程 拟 合 
度 均 达到 了 极 显著 水 平 。 其 中 泡 泡 刺 (Witraria 
sphaerocarpa ) 、 旱 (Artemisia xerophytica) 、 梭 梭 和 和 
树 柳 的 最 佳 预 测 变 量 为 V; 沙 拐 束 、 膜 果 有 麻黄 PRR 
木 红 砂 、 合 头 葵 和 霸王 的 最 佳 预 测 变 量 为 4, 其 最 
优异 速生 长 方程 见 表 1。 

根据 每 个 样 方 调查 数据 和 灌木 异 速生 长 方程 
以 及 小 样 方 草本 植物 地 上 生物 量 调查 数据 ,对 每 个 
样 方 的 地 上 生物 量 进行 计算 ,计算 结果 用 于 荒漠 植 
被 地 上 生物 量 遥 感 估算 。 


表 1 灌木 植物 的 最 优异 速生 长 方程 


Tab.1 Bestallometric model of shrub plants 


物种 预测 变量 异 速生 长 方程 R MAE/g 植株 个 数 

MMI (Nitraria sphaerocarpa) V y2131.27 y**! 0.87 29.16 n=15(P<0.0001 ) 
IUE (Artemisia xerophytica ) V y-393.99 V°” 0.86 87.54 n-14(P«0.0001 ) 
梭 梭 (Haloxylon ammodendron ) V y=688.38V°" 0.78 235.25 n=15(P<0.0001 ) 
FEM (Tamarix ramosissima ) y y=533.09V°" 0.73 355.42 n=15(P<0.0001 ) 
WFR ( Calligonum mongolicum) A y=18.70A-2.71 0.99 2.89 n=15( P<0.005) 
膜 果 麻黄 (Ephedra przewalskii) A y=21.49A-12.49 0.97 8.42 n=13(P<0.005) 
TUR ZR (Gymnocarpos przewalskii ) A y712.844-2.36 0.99 2.07 n=9( P<0.005) 
红 砂 (Reaumuria soongarica ) A y=11.76A-1.10 0.74 5.34 n=23(P<0.005 ) 
EKHE ( Sympegma regelii ) A y77.944-1.58 0.83 4.80 n=28( P<0.005) 
霸王 (Zygophyllum xanthoxylon ) A y=7.93A-1.84 0.96 1.25 n=10(P<0.005 ) 


注 : 尼 为 决定 系数 ;MAE 为 平均 绝对 误差 ; 为 植株 体积 (cm ) ;4 为 冠 面积 (cm ) ; y 为 植被 地 上 生物 量 佑 算 值 (g) ;7 为 样本 个 数 ; 忆 为 显著 性 


水 平 。 
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2.3 遥感 数据 处 理 

2.3.1 Landsat-OLI i$ £i 4 4& 4E 3€. Landsat-OLI 数 
据 由 美国 地 质 勘 探 局 网 站 下 载 (https://earthexplorer. 
usgs.gov/) ,获取 时 间 是 2017 年 9 月 26 日 , 履 盖 于 阿 
尔 金山 北 芒 山 前 戈壁 区 ,空间 分 辩 率 30 m。Land- 
sat-OLI 数 据 包 含 可 见 光 (VIS)、 近 红外 (NIR) 和 短 红 
外 (SWIR) 等 9 个 光谱 波段 。Landsat-0LI 数 据 是 经 
过 陆地 卫星 生态 系统 干扰 自 适应 处 理 系 统 处 理 后 
的 地 表 反 射 率 产品 ,可 直接 用 于 遥感 特征 参数 的 计 
算 。 

2.3.2 WorldView-3 i£ £i 4& 4& ^E 3€. WorldView-3 
BC 48 TK CAN [8] Fy 2016 Æ 6 H 28 H , TEE Thi FR 381 
km’, Z [8] 4) EAS 2.5 m, 包含 蓝光 (Blue) , 2856 
(Green) 红 光 (RED) 和 NIR 波 段 。 根 据 课题 组 在 地 
面 样 方 附近 山药 地 带 设 置 的 自动 气象 站 ( 距 World- 
View-3 遥感 影像 覆盖 区 域 西部 边缘 约 10 km) 的 观 
测 数据 ,2016 年 降水 量 75.5 mm,2017 年 降水 量 56.8 
mm ; 该 地 区 气候 极端 干旱 , 苇 漠 植物 生长 缓慢 ,1a 
内 植被 生物 量变 化 很 小 。 遥 感 影像 数据 获取 与 地 
面 样 方 调查 之 间 的 时 间 差 异 不 会 给 植被 地 上 生物 
量 估算 结果 带 来 显著 影响 。 

在 ENVI 5.3 环境 中 ,利用 数据 元 文件 对 原始 的 

WorldView-3 数据 进行 辐射 定 标 ; 通 过 高 光谱 遥感 
数据 的 大 气 纠正 模块 (FLAASH) 对 WorldView-3 数 
据 进 行 大 气 校正 ;利用 影像 的 正 射 校 正文 件 (PRC ) 
和 全 球 定 位 系统 (GPS ) 差 分 定位 对 影像 进行 几何 
精 校正 ,误差 可 以 控制 在 2m 以 内 ;提取 遥感 特征 
BR, 
2.3.3 遥感 特征 参数 计算 ”选取 单反 射 率 波段 和 植 
被 指数 作为 遥感 特征 参数 ,用 于 荒漠 植被 覆盖 区 的 
提取 和 植被 地 上 生物 量 回归 模型 的 构建 。 单 反射 
率 波 段 中 RED 和 NIR 波段 对 荒漠 植被 的 敏感 性 较 
强 呈 ;植被 指数 中 RVI.NDVI 和 MSAVI 对 荒漠 植 
被 的 敏感 性 更 强 2 5 。 因 此 ,本 研究 采用 的 遥感 特 
征 参 数 有 5 个 , 即 RED NIR ,RVI,NDVI fll MSAVI, 

公式 (3) 公式 (4) 和 公式 (5) 为 植被 指数 的 计 
算 公 式 , 式 中 RED 和 NIR 分 别 表示 红 光 和 近 红 光 反 
射 率 波段 。 


RVI= NIR/RED (3) 
NDVI=(NIR - RED)/(NIR + RED) (4) 


MSAVI= [2NIR +1-_/QNIR + If - 8(NIR - RED) |/2 (5) 


2.3.4 A Aik BH dg eZ 202049 H 26—30 
日 ,采用 DJI 精 录 4(RTK) 无 人 机 在 WorldView-3 数 
据 覆 盖 区 内 设置 了 一 条 样 带 进行 详细 调查 ,覆盖 区 
域内 的 荒漠 植被 具有 代表 性 。 该 样 带 南 北向 长 11 
km ,东西 向 宽 4km, 由 44 个 1 kmxl km 的 飞行 单元 
组 成 (图 1)。 为 避免 太阳 入 射 角度 太 小 而 导致 的 地 
物 阴影 拖 中 现象 ,拍摄 时 间 选 取 在 北京 时 间 09:00 一 
16:00 之 间 。 飞 行 作业 期 间 天 气 晴 好 , 云 量 极 低 , 风 
速 很 小 ,满足 无 人 机 摄影 测量 平台 工作 的 要 求 ; 飞 
行 高 度 为 200 m, 飞 行 速度 9 m-s”, ALIA] FER 8096, 
35 [n] TE € 75% ;空间 分 辩 率 0.05 m。 数 据 的 平面 
坐标 系统 为 WGS84 坐标 系 ,通用 墨 卡尔 投影 。 无 人 
机 影像 覆盖 区 域 的 荒漠 植被 具有 代表 性 ,灌木 冠 幅 
一 般 大 于 20 cm ,在 无 人 机 影像 上 可 被 精准 识别 。 

2.4 植被 覆盖 区 的 提取 与 精度 评价 

2.4.1 植被 覆盖 区 提取 原理 “利用 像 元 二 分 法 提取 
植被 覆盖 区 ,以 提高 植被 地 上 生物 量 模型 估算 精 
度 'W。 其 原理 是 假设 地 表 信 息 由 纯 植 被 和 纯 土 壤 
言 息 构 成 , 则 与 之 相对 应 的 像 元 光谱 信息 也 由 一 定 
履 盖 比例 的 纯 植 被 像 元 与 纯 土 壤 像 元 构成 。 表 
示 为 : 


S=FVCxS y +l- FVO X Si (6) 
XAR (6 ) HEFT AE He, FEE AEE ii EIAS 
FVC= (Sa SMS x (7) 


式 中 :FVC 为 像 元 中 纯 植 被 覆盖 面积 的 比例 (植被 
T5 s RE) S 为 像 元 的 遥感 特征 参数 值 ;9 和 SS 分别 
为 像 元 中 纯 植 被 和 纯 土 壤 像 元 的 遥感 特征 参数 值 。 
2.4.2 植被 和 土壤 分 类 ”利用 目 视 解 译 从 无 人 机 遥 
感 影像 中 提取 纯 土壤 像 元 分 布 区 ,其 他 像 元 对 应 的 
区 域 则 为 植被 覆盖 区 ”“”。 通 过 非 监 督 分 类 方法 对 
所 有 无 人 机 影像 进行 植被 和 土壤 分 类 ,拼接 分 类 结 
果 得 到 植被 和 土壤 分 类 图 。 图 2 为 一 个 地 面 调查 样 
方 的 无 人 机 影像 及 其 分 类 结果 , 样 方 中 心 地 理 坐标 
为 39"20'34"N .91°59'4"E, KA 30 mx30 m。 从 图 
2 可 以 看 出 ,无 人 机 影像 能 够 清晰 地 反映 样 方 内 低 
矮 稀 玻 荒 漠 植 被 的 冠 层 信息 与 纯 土 壤 像 元 信息 。 
Landsat-OLI 和 WorldView-3 影 像 的 植被 覆盖 区 
的 提取 方法 一 致 。 采 用 非 监 督 分 类 方法 分 别 对 遥 
感 影像 进行 植被 和 土壤 分 类 ; 掩 膜 提 取 土 壤 覆 盖 
区 ;统计 土壤 覆盖 区 遥感 特征 参数 的 累计 频率 。 对 
照 无 人 机 影像 数据 的 植被 和 土壤 分 类 结果 并 参考 
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TA 


图 例 

mm 植被 覆盖 区 
O 土壤 覆盖 区 
0 0.006 km 
L 


TE PED EN 17% 
图 2 样 方 无 人 机 影像 及 其 分 类 结果 
Fig.2 Plot UAV image and its classification results 


相关 文献 ,经 过 反复 调试 ,选取 频率 累计 表 中 累计 
频率 为 5% 的 遥感 特征 参数 作为 纯 土 壤 像 元 Sm 
值 ,遥感 特征 参数 累计 频率 >5% 的 像 元 为 植被 覆 
Km! , 

2.4.3 植被 覆盖 区 提取 与 验证 ”以 无 人 机 遥感 数据 
的 植被 分 类 结果 作为 实测 值 ,以 遥感 影像 的 分 类 结 
果 作 为 预测 值 ,利用 Kappa 系数 对 分 类 结果 进行 评 
价 , 计 算 公 式 如 下 : 


Kappa= 


P,-P. 
T-P. (8) 
式 中 :Po 为 实测 值 与 预测 值 相同 的 像 元 数量 ;P. 为 实 
测 值 与 预测 值 不 同 的 像 元 数量 。 

分 别 统计 Landsat-OLI FN WorldView-3 遥感 数 
据 的 5 个 特征 参数 影像 的 土壤 覆盖 区 像 元 信息 。 土 
壤 覆 盖 区 内 累计 频率 为 5% 的 特征 参数 值 代表 纯 土 
坏 像 元 信息 ,Landsat-OLI 遥 感 数据 的 5 个 特征 参数 
对 应 的 Su 值 分 别 为 0.27.、0.31、1.13 0.06 和 0.05， 
WorldView-3 遥感 数据 的 5 个 特征 参数 对 应 的 Sm 值 
分 别 为 0.26.0.22 .1.11 .0.04 和 0.03。 累 计 频 率 >5% 
的 像 元 覆盖 区 域 为 植被 覆盖 区 。 利 用 无 人 机 还 感 
数据 分 别 对 Landsat-OLI[ 和 WorldView-3 遥感 数据 
的 植被 和 土壤 分 类 结果 进行 精度 评价 。 由 表 2 可 
知 ,Landsat-OLI 的 RVI 分 类 效果 最 优 ,其 Kappa 系数 


H 0.77; WorldView-3 的 NDVI 分 类 效果 最 优 , 其 
Kappa 系数 为 0.95。 分 别 采用 最 优 植被 指数 RVI 和 
NDVI 对 研究 区 进行 植被 覆盖 区 的 提取 ,World- 
View- 3 影像 覆盖 区 域 的 植被 覆盖 区 面积 为 223 
km' ,Landsat-OLI 影 像 覆 盖 区 域 的 植被 覆盖 区 面积 
H 3592 km?。 在 无 人 机 影像 覆盖 区 随机 选取 200 个 
像 元 对 植被 覆盖 区 提取 结果 进行 验证 ,Landsat-OLI 
影像 的 植被 覆盖 区 提取 精度 为 84.54% , WorldView- 
3 影像 的 植被 覆盖 区 提取 精度 为 91.26% 。 
2.5 植被 地 上 生物 量 遥 感 模型 的 构建 与 验证 

将 研究 区 内 100 个 样 方 的 中 心 点 GPS 坐标 羞 加 
于 Landsat-OLI 影 像 上 ,从 中 随机 选取 70%(n=70) 的 
样 方 用 于 模型 构建 ,其 余 30%(n=30) 的 样 方 用 于 模 
型 验证 。 建 立 70 个 样 方 的 植被 地 上 生物 量 数据 与 
Landsat-OLI 遥 感 影像 上 对 应 像 元 的 光谱 特征 值 之 
间 的 回归 模型 。100 个 样 方 中 有 21 个 样 方 分 布 于 
WorldView-3 影像 覆盖 区 ,将 这 21 个 样 方 均 分 为 84 
个 15 mx15 m 的 小 样 方 。 样 方 内 所 有 植株 的 相对 位 
置 已 知 ,可 以 准确 计算 每 个 小 样 方 的 植被 地 上 生物 
量 , 将 其 中 心 点 的 GPS 坐标 车 加 于 WorldView-3 影 
像 上 ,随机 选取 70%(n=58) 的 小 样 方 ,建立 小 样 方 
地 上 生物 量 数 据 与 WorldView-3 遥感 特征 参数 之 间 
的 回归 模型 ,剩余 30% 的 小 样 方 (n=26) 用 于 回归 模 


表 2 植被 和 土壤 分 类 植被 指数 


Tab.2 Vegetation indices for vegetation and soil classification 


m MUR Landsat-OLI WorldView-3 
遥感 特征 参数 
RED NIR RVI NDVI MSAVI RED NIR RVI NDVI MSAVI 
Kappa 系数 0.72 0.72 0.77 0.70 0.75 0.82 0.85 0.90 0.95 0.93 


TE: RED 为 红 光 反射 率 波段 ;NIR 为 近 红 外 反射 率 波段 ;RVI 为 比值 植被 指数 ;NDVI 为 归 一 化 植被 指数 ; MSAVI 为 修正 土壤 调节 植被 指数 。 
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叶 静 又 等 极 干旱 区 稀疏 荒漠 植被 地 上 生物 量 通 感 估算 


型 验证 。Townshend 等 研究 表明 ,遥感 数据 验证 
ERED AV) ESF (0482) ,其 中 巨 为 样 方 大 小 
(m) ,下 为 像 元 大 小 Am) , 工 为 像 元 几何 精度 (m)。 
WorldView-3 遥感 数据 的 像 元 大 小 为 2.5m ,几何 精 
度 2mm 以 内 , 则 最 佳 样 方 大 小 为 12.5m, 其 与 15 mx 
15 m 的 小 样 方 大 小 相近 。 因 此 , 样 方 数据 与 遥感 像 
元 之 间 存 在 匹配 关系 。 

SLR 、LASSO、RR 模 型 常 被 应 用 于 植被 地 上 生 
物 量 的 预测 玉 。 在 R 语 言 环境 中 构建 回归 模型 ， 
利用 “glmnet" 包 构建 遥感 特征 参数 与 植被 地 上 生物 
量 之 间 的 回归 模型 ,并 利用 “sperrorest” 包 进行 模型 
的 交叉 验证 "。 交 又 验 证 数据 集 被 随机 分 配 至 G 
(G=10) 组 ,交叉 验证 重复 10 次 ,选取 均 方 根 误 差 
(RMSE) .标准 偏差 (SD) 和 均 方 根 误差 百分比 (RM- 
SErel) 作 为 模型 精度 的 指标 ,用 于 对 比分 析 模 型 的 回 
归 效 果 。RMSE 用 来 衡量 预测 值 与 实测 值 之 间 的 偏 
差 ;SD 为 各 数据 偏离 平均 数 的 距离 ;RMSErel 为 均 方 
根 的 偏 移 。3 项 指标 越 小 表示 模型 的 误差 越 小 。 


(9) 
(10) 
RMSErel = ET x 100% (11) 


式 中 : y, HP i MEET LE ES T Cg m 7) ; 
7 为 第 i 个 样 方 植被 地 上 生物 量 预测 值 (gm?); y 
为 实测 植被 地 上 生物 量 的 平均 值 (gm 7) ;nn 为 模型 
中 使 用 的 样 方 个 数 。 
2.6 植被 地 上 生物 量 遥 感 估 算 结 果 的 修正 

构建 WorldView-3 和 Landsat-OLI 影 像 植 被 地 
上 生物 量 分 布 数据 之 间 的 线性 关系 ,通过 该 线性 关 
系 对 基于 Landsat-0LI 影 像 得 到 的 阿尔 金山 北 芳 山 
前 戈壁 区 的 车 漠 植 被 地 上 生物 量 分 布 图 进行 修正 。 


3 结果 与 分 析 


3.1 估算 结果 验证 

采用 交叉 验证 方法 对 所 有 样 方 的 估算 结果 进 
行 对 比分 析 ( 表 3)。WorldView-3 数 据 的 反 演 效果 
明显 优 于 Landsat-OLI, RMSE 平 均 降 低 58.65 g «m^? , 
RMSErel 平 均 降 低 了 72.27%。LASSO .RR 和 SLR 模 


表 3 植被 地 上 生物 量 预测 结果 交叉 验证 
Tab.3 Result of cross-validated aboveground 


biomass in the samples 


模型 类 型 RMSE/g:m? SD/g:m* | RMSErel/% 
Landsat-OLI LASSO 95.36 42.84 115.81 
WorldView-3 LASSO 32.99 10.27 45.26 
Landsat-OLI RR. 94.15 40.99 114.34 
WorldView-3 RR 33.42 9.81 44.68 
Landsat-OLI SLR 95.92 43.39 134.94 
WorldView-3 SLR 43.07 21.93 58.33 


iE : Landsat- OLI 为 Landsat- OLI f£ /& $$ ; WorldView- 3 7j World- 
View-3 fidi ; LASSO 为 套 索 回归 ;RR 为 岭 回归 ;SLR 为 逐步 线性 
"IH ; RMSE 为 均 方 根 误 差 ;SD 为 标准 偏差 ;RMSErel 为 均 方 根 误 
差 百 分 比 。 


型 中 ,RR 模型 的 估算 效果 最 好 ,其 中 WorldView-3 
RR 模型 的 RMSErel 为 44.68% , 略 优 于 WorldView-3 
SLR 模型 的 58.33% 和 WorldView-3 LASSO 模型 的 
45.26% ; Landsat-OLI RR 模型 的 RMSFrel 7j 114.34% , 
略 优 于 Landsat-OLI SLR 模型 的 134.94% 和 Landsat- 
OLI LASSO 模 型 的 115.81%。 

3.2 估算 结果 分 析 

利用 WorldView-3 RR 和 Landsat-OLI RR 模型 
估算 植被 地 上 生物 量 , 得 到 模型 预测 值 与 实测 值 的 
1:1 线性 图 (图 3)。 利 用 WorldView-3 RR 模型 反 演 
得 到 的 植被 地 上 生物 量 预测 值 与 样 方 实测 值 更 巾 
jE F 1:1 2B , SEF WorldView-3 数据 的 估算 效果 明显 
=F Landsat-OLI RR. Landsat-OLI RR 模型 植被 地 
上 生物 量 估算 值 基 本 上 不 随 实测 值 的 变化 而 发 生 
ARAL, , REH Landsat-OLI Egg A He MT T fibt oc D A BK 
的 敏感 性 非常 低 ,在 实际 应 用 中 会 使 植被 地 上 生物 
量 估算 结果 出 现 较 大 偏差 。 

在 WorldView-3 数 据 覆 盖 区 内 随机 选取 500 个 
像 元 ,每 个 像 元 之 间 的 距离 不 小 于 100 m, 构建 
WorldView-3 RR 和 Landsat-OLI RR 植被 地 上 生物 
量 估算 数据 之 间 的 线性 关系 ,对 Landsat-OLI RR fH 
被 地 上 生物 量 估算 结果 进行 修正 。WorldView-3 RR 
5j Landsat-OLI RR 估算 值 之 间 的 关系 式 为 (图 4): 

OLI= 0.1281WV +71.964 (12) 
XP OLI 5; WV 4} 5!) W Landsat- OLI RR 5j World- 
View-3 RR 模型 的 植被 地 上 生物 量 预测 值 (g.m”)。 
3.3 植被 地 上 生物 量 遥 感 估算 

利用 WorldView-3 RR 估算 WorldView-3 影像 
履 盖 区 的 植被 地 上 生物 量 ,其 均值 为 122.43 gm? 
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500 r (a) WorldView-3 RR 模型 
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(b) Landsat-OLIRR 模 型 


预测 值 /gm 


实测 值 /g-m? 


图 3 植被 地 上 生物 量 实测 值 和 预测 值 的 1:1 线 性 图 


Fig.3 Scatterplots of measured and predicted aboveground biomass 
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注 :0LI 与 WV 分 别 为 Landsat-OLI RR 5j WorldView-3 RR 模型 的 植 
被 地 上 生物 量 预测 值 ; 尼 为 决定 系数 。 
图 4 植被 地 上 生物 量 回归 分 析 


Fig. 4 Regression analysis of aboveground biomass 


(a)WorldView-3 RR 模型 估算 结果 


(c) 修正 后 的 Landsat-OLI RR 模型 估算 结果 


图 5 WorldView—3 遥感 影像 覆盖 区 域 植 被 地 上 生物 量 佑 算 结 果 


(RMSErel=44.68% ) ( d 5a) ;利用 Landsat-OLI RR fi; 
fi WorldView-3 影像 覆盖 区 的 植被 地 上 生物 量 ,其 
均值 为 84.59 g - m? (RMSErel-114.3496 ) (图 5b) 。 
Landsat-OLI RR 的 RMSErel Lt WorldView-3 RR 高 出 
69.66% 。 利 用 公式 (12) 对 Landsat-OLI RR 模型 估 
算 结果 修正 ,修正 的 植被 地 上 生物 量 均值 为 90.35 
g*m "(RMSErel=83.16%)( 图 5c)。 修 正方 程 在 一 定 
程度 上 提高 了 Landsat-OLI RR 模型 的 估算 精度 ， 
RMSErel 降 低 了 31.18%。 对 比 图 5b 和 图 5c 可见, 修 
正 后 的 植被 地 上 生物 量 分 布 更 加 “精细 ”。 洪 积 扇 
戈壁 表面 因 洪 水 而 沟 密 纵横 , 受 微 地 形 影 响 植 被 分 
布 空间 异 质 性 较 高 ,图 5c 对 植被 地 上 生物 量 分 布 的 


(b) Landsat-OLI RR 模型 估算 结果 
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Fig. 5 Estimation results of aboveground biomass in the area covered by WorldView-3 remote sensing data 
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(a) Landsat-OLI RR 模型 估算 结果 
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6 植被 地 上 生物 量 估算 结果 


Fig. 6 Estimation results of aboveground biomass 


空间 异 质 性 刻画 得 更 加 清晰 。 

利用 Landsat-OLI RR 模型 估算 研究 区 植被 地 上 
生物 量 ,其 均值 为 85.28 g +m ?(RMSErel=84.93% ) 
(图 6a)。 修 正 后 的 植被 生物 量 均值 为 118.28 gm? 
(RMSErel=71.51%)( 图 6b)。 从 植被 地 上 生物 量 分 
布 来 看 ,图 6b 的 植被 地 上 生物 量 分 布 格局 与 地 面 调 
查 数据 更 为 符合 ,南部 海拔 较 高 的 山药 地 区 植被 地 
上 生物 量 较 高 ,而 北部 接近 于 沙漠 边缘 的 海拔 较 低 
区 域 植被 地 上 生物 量 较 低 。 对 比 两 幅 地 上 生物 
量 分 布 图 可 以 发 现 , 在 植被 地 上 生物 量 较 低 的 
93°00'E 附 近 地 区 和 植被 地 上 生物 量 普 遍 较 高 的 
94°20' 忆 附近 地 区 ,修正 前 后 分 别 出 现 了 高 值 变 低 
和 低 值 变 高 的 现象 ,表明 低 植被 地 上 生物 量 区 域 被 
高 估 和 高 植被 地 上 生物 量 区 域 被 低估 的 问题 得 到 
了 一 定 程度 的 修正 。 
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C1) RENCE AS BA BER D UU t Pr HL 
地 上 生物 量 回归 模型 估算 精度 的 主要 因素 。 与 中 
空间 分 辨 率 Landsat-OLI 遥 感 数据 相 比 , 基 于 高 空间 
AT BER World View—3 遥感 数据 建立 的 稀 玻 荒漠 植被 
地 上 生物 量 估算 模型 估算 精度 更 高 。 

(2) 传统 的 植被 地 上 生物 量 回归 模型 多 是 利用 
植被 指数 等 遥感 特征 参数 与 地 面 数 据 的 相关 关系 
直接 构建 。 虽 然 采 用 该 方法 可 能 会 得 到 较 高 的 估 
算 精 度 ,但 是 遥感 特征 参数 之 间 的 共 线 问题 会 使 模 


型 的 稳定 性 降低 。 采 用 非 线性 的 多 元 回归 方法 
SLR ,LASSO 和 RR 构 建 遥感 特征 参数 与 地 面 样 方 数 
据 之 间 的 回归 模型 可 以 有 效 提高 模型 的 稳定 性 。 

(3) 提取 稀 玻 荒漠 植被 信息 时 ,采用 高 空间 分 
辩 率 遥感 数据 的 佑 算 结 果 对 中 空间 分 辩 率 遥感 数 
据 的 估算 结果 进行 修正 ,可 以 有 效 提高 中 空间 分 辨 
率 遥 感 数据 对 稀 玻 荒漠 植被 信息 的 提取 精度 ,并 可 
以 在 具有 相同 自然 条 件 的 较 大 尺度 上 应 用 。 
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Estimation of aboveground biomass of sparse desert vegetation 


based on remote sensing techniques in hyper-arid area 
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Abstract: The aboveground biomass (AGB) is an essential indicator of desert ecosystem health evaluation and 
desertification monitoring in arid areas. The estimation accuracy of the general remote sensing model is reduced 
because of the influence of small and sparse desert vegetation, low AGB per unit area, and background soil 
information. Improving biomass estimation accuracy of sparse desert vegetation can provide essential parameters 
for terrestrial ecosystem carbon storage estimation and carbon cycle research. It can also provide technical 
support for desertification monitoring. Herein, the Gobi area in front of the Altun Mountains in the south of 
Kumtag Desert, northwest China with an extremely arid climate is considered as the research area. The vegetation 
coverage area was accurately extracted using UAV image data. The estimation results based on the high-spatial 
resolution WorldView-3 remote sensing data were used to correct the estimation results based on the medium 
spatial resolution Landsat- OLI remote sensing data. Then, stepwise linear regression (SLR), LASSO regression, 
and ridge regression (RR) models were used to estimate the AGB of sparse desert vegetation. The results showed 
that: (1) High-spatial resolution satellite remote sensing data can accurately describe the spatial and temporal 
distribution characteristics of vegetation in arid areas. The RMSErel mean value of WorldView-3 remote sensing 
data was 44.68%. (2) The RR model had a better regression effect than SLR and LASSO models. The RMSErel 
of the WorldView-3 RR model was 44.68%, better than the SLR model (58.33%) and LASSO model (45.26%). 
Additionally, the RMSErel of the Landsat- OLI RR model was 114.3495, better than the SLR model (134.9494) 
and LASSO model (115.8196). (3) Using the WorldView-3 remote sensing data as an intermediate medium 
between ground and Landsat-OLI remote sensing data, the AGB estimation results of desert vegetation based on 
Landsat- OLI remote sensing data were corrected, and the RMSErel decreased by 31.18%. The spatial resolution 
of remote sensing data is the main factor affecting the estimation accuracy of the sparse desert vegetation AGB 
regression model. The main reason for the decrease in the estimation accuracy of AGB is the decrease in the 
sensitivity of Landsat- OLI remote sensing data to sparse desert vegetation information. The nonlinear multiple 
regression methods, SLR, LASSO, and RR, are used to build the regression model between remote sensing 
characteristic parameters and ground sample data, effectively improving model stability. When extracting sparse 
desert vegetation information, using high-spatial resolution remote sensing data as a medium between ground and 
medium spatial resolution remote sensing data can effectively improve the extraction accuracy of sparse desert 
vegetation information from medium spatial resolution remote sensing data. 


Key words: remote sensing eigenvalue; sparse desert vegetation; aboveground biomass; regression model 


